The Resonant Exchange Qubit 
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We introduce a solid-state qubit in which exchange interactions among confined electrons provide 
both the static longitudinal field and the oscillatory transverse field, allowing rapid and full qubit 
control via rf gate- voltage pulses. We demonstrate two- axis control at a detuning sweet-spot, where 
leakage due to hyperfine coupling is suppressed by the large exchange gap. A 7r/2-gate time of 2.5 ns 
and a coherence time of 19 /is, using multi-pulse echo, are also demonstrated. Model calculations 
that include effects of hyperfine noise are in excellent quantitative agreement with experiment. 

PACS numbers: 



As originally conceived, the two-level system that 
forms the basis of the semiconductor spin qubit is the 
electron spin itself, with pulsed exchange between two 
confined electrons forming a two-qubit gate [1 . Gener- 
alizations to two-electron [2H6] and three-electron |7HT3] 
qubits make use of multi-electron states as the quantum 
two-level system. These qubits offer ease of initialization, 
control, and readout, or speed of operation, in exchange 
for the complexity of controlling more than one electron 
per qubit. An attractive feature of the original single- 
spin proposal is that qubit rotations are implemented 
as Rabi processes, driven by a small resonant transverse 
field, rather than Larmor processes, which use pulsed 
Larmor precession around larger nonpar allel fields. Rabi 
rotations allow narrow-band wiring away from dc, preces- 
sion rates controlled by the amplitude of the oscillatory 
field, and straightforward two- axis control (needed for ar- 
bitrary transformations) implemented using the phase of 
the oscillatory field [lHHS]. 

In this Letter, we introduce a new quantum-dot-based 
qubit — the resonant exchange qubit — that captures the 
best features of previous incarnations, with qubit rota- 
tions via Rabi nutation using gate-controlled exchange 
both for the static longitudinal field and the oscillatory 
transverse field, as described in Ref. [16]. The large ex- 
change field suppresses leakage from the qubit space. 
However, because rotations are driven by a resonant 
transverse field, the large longitudinal field does not im- 
pose unrealistically fast evolution between qubit states. 
Moreover, the qubit is operated at a "sweet spot" of the 
exchange gap, making it insensitive to first order to elec- 
trical noise in the detuning parameter p!6UT9] . 

The resonant exchange qubit was realized in a triple 
quantum dot formed by surface gates 110 nm above a 
two-dimensional electron gas (density 2.6 x 10^^ m~^, 
mobility 43 m^/Vs) in a GaAs/Alo.sGao.rAs heterostruc- 
ture [see Fig. [ija)]. Gate voltages Vi and Vr controlled 
detuning, e = {Vr - V^)/2 - {Vi - V^°)/2, measured rel- 
ative to the center of the 111 charge region, while Vm 



controlled the size of the 111 region (111 and other num- 
ber triplets denote the charge occupancy of the triple 
dot) [20 . An adjacent multi-electron quantum dot oper- 
ated in Coulomb blockade regime served as a radio fre- 
quency (rf ) charge sensor [HI [22] . 

Tunneling between adjacent quantum dots gives two 
exchange splittings, Ji{s)^ associated with the electron 
pair in the left and middle dots, and Jr{^)^ associated 
with the electron pair in middle and right dots. Away 
from zero detuning, defined as the center of 111, the 
qubit ground state, |0) = ^(|m) + lltt) - 2|nt)), 
connects continuously to a singlet state of the left pair, 
\Si) = ^(Itlt) - lltt)) in charge state 201, and to a 
singlet state of the right pair, l^^) = T^dttl) - Itit)) 
in charge state 102. [see Fig. 2(a)]. The excited qubit 
state, |1) = -^(im) - lltt)), maps into triplet states 
that, in contrast to the singlets, cannot tunnel into charge 
states 201 or 102. This allows the qubit state to be de- 




FIG. 1: (a) False color micrograph of lithographically iden- 
tical device with dot locations depicted; gates are marked in 
yellow. Gate voltages, Vi and Vr, set the charge occupancy of 
left and right dot as well as the detuning, e of the qubit. A 
neighboring sensor quantum dot is indicated with a larger cir- 
cle, (b) Triple dot charge occupancy N^ Nm Nr as a function 
of Vi and Vr in and near the 111 regime; e — {Vr — Vr)/2 — 
{Vi - V^°)/2, 5=(Vr- V^)/2 + {Vi- V^)I2 + ^{V^ - V^). 
Measurements give 7^3. The operating position is marked 
with a star, which is larger than the amplitude of voltage 
fluctuations used in rotations. 




latory voltage to gate Vi, which moves the operating point 
around e = 0, in turn creating an oscihatory transverse 
field Jx [see Fig.[2|b)]. When the oscillation frequency lo 
matches the longitudinal exchange frequency, Jz/h [see 
Fig.j2Vb)], the qubit nutates between |0) and |1). Fig- 
ure ^c) maps the positions of the |Q+) anti-crossings 
with the lower qubit branch as a function of field and de- 
tuning without applied microwaves, along with a model 
calculation of the exchange splittings Ji and J^. This 
spectroscopy is performed by preparing a l^^) state in 
102, then pulsing into 111 for 300 ns before returning to 
102 to project the resulting state back onto l^^). 

The data in Fig. [2|d) shows two features, a vertical 
line corresponding to the crossing of \Q-\-) and the center 
of the lower qubit branch (circle), and a curved feature 
refiecting a driven oscillation between qubit states |0) and 
|1), marked with a star. The curved feature shows that 
the qubit splitting is controlled by gate voltage Vm^ here 
covering a range from 200 MHz to 2 GHz. Using fast 
gating, we have demonstrated control of this frequency 
on nanosecond time scales. The dashed line in Fig. |2|d) 
is a model of uj{Vm) that assumes a linear dependence of 
Ji and Jr on Vm- 

The resonant exchange qubit can be modeled by the 
Hamiltonian, 



FIG. 2: (a) Energy level diagram for a constant 6. Charge 
transitions are marked with circles (qubit- |Q+)), triangles 
(qubit- 1 Q+), photon), squares (qubit- |Q+), photon), and a 
star (|0)-|1)) transitions, (b) Schematic of the effects of Ji and 
Jr on the qubit Bloch sphere (c) The qubit- 1 (5+) anti-crossing 
is mapped out in magnetic field and detuning without an ex- 
citation. Dashed line is a model of the exchange splitting for 
equal tunnel couplings, (d) A sweep of the middle plunger 
gate at s = mV and fixed field of 575 mT, demonstrating 
control of the main qubit transition. The dashed curve is a 
model of Ji + Jr as a function of Sq — {Vm — V^)/2 [23] . 

(e) At a fixed field of 310 mT, detuning and microwave burst 
frequency are swept to trace out the spectroscopy of the qubit. 

(f ) A model of qubit evolution in the presence of a microwave 
excitation and magnetic field gradients between dots in the 
longitudinal and transverse directions. 



tected with a charge sensor the distinguishes 201, 111, 
and 102. A third state, |(3+) = |ttt), intersects the 
qubit ground state at two anti-crossings whose posi- 
tion depends on Zeeman splitting from an external mag- 
netic field. By sweeping the magnetic field, the qubit 
ground-state energy can be measured as a function of 
detuning [Figs. |2|a,c)]. The fourth state in Fig. |2|a), 
\Q) = Tjdtti) + Itit) + litt)), is separated from the 
qubit states by a sizable gap (half the separation between 
|0) and |1)), suppressing leakage out of the qubit space. 
The gap to \Q) is deliberately kept large by setting tun- 
neling rates, hence Ji and J^, to be large throughout the 
111 charge region. 

Qubit rotations are implemented by applying an oscil- 



n{e) = -JzCJz/2 - JxCTx/2, 



(1) 



^{Me)-Ji{e)), 



where Jz = ^{Ji{e) -\- Jr{e)) and Jx — 2 
where a^ = |0)(0| - |1)(1| and ax = |0)(1| + |1)(0| are 
the Pauli operators of the qubit [see Fig. [2 |b)]. Ex- 
change fields Ji{e) = -(e + eo)/2 + y/t^ ^ {e -^ e^)^ /A 
and Jrie) = {e — Sq)/'! + ^/f^ + (s — s^Y /A are mod- 
eled in terms of the tunnel coupling, t, which is taken 
to be the same for both the 201-111 and 111-102 tran- 
sitions, and i^O' ^^^ detunings of these charge transi- 
tions. At £ = 0, this gives dJ^/de = and dJx/de = 

^(1 — SqI ^U'^ +^0^)- ^^^ smah detuning, e <^ e^, Jz 
is unchanged to first order while Jx ^ e. This system 
is equivalent to a spin- 1/2 in a large static field with a 
small transverse field. While Jz is insensitive to detuning 
noise to first order, it is not insensitive to noise on gate 
Vm or other gates. However, other gates, including Kn, 
do not need to operate at high frequency, and so can be 
heavily filtered. 

In Fig. |3J \Sr) is prepared in 102 and adiabatically 
evolved to |0) at e = 0, taking care to move rapidly 
through the \Q-\-) anti-crossing. A microwave burst is 
then applied to Vi for a time tb before returning adia- 
batically to 102 for measurement. The color plot shows 
the probability, Pq^ of detecting the ground state through 
a charge measurement (see Sec. I, Supplemental Mate- 
rial.) By sweeping frequency and power, we see patterns 
characteristic of Rabi nutations subject to low frequency 
noise in the splitting frequency, cjqi due to hyperfine 
gradients (see Sec. VI, Supplemental Information). In 
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FIG. 3: (a) A Rabi nutation for a -45 dBm (0.45 mV) 
excitation on the left plunger gate with the detuning biased 
to the center of the transition [(•) in Fig. [2[a,e)] for a time 
tb. (b) A Rabi nutation with a 0.355 GHz excitation on the 
left plunger gate with the detuning biased to the center of 
the transition, (insets) A model of this nutation using the 
exchange profile from Fig. l2{f ) and fluctuating longitudinal 
magnetic field gradients. 
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FIG. 4: (a) A schematic of the detuning during a two-pulse 
sequence, where the first pulse is an X rotation and the second 
pulse is a rotation around an angle set by the relative phase of 
the carrier, $, as depicted on the Bloch sphere, (b) The qubit 
readout for a rotation about X, followed by a ^ rotation 
about an axis $, for three different $'s. The solid lines are 
fits to the model in Fig. [2|c,d,f ) and the insets of Fig. Is] 



the rotating frame, the amplitude of the oscillation gives 
the strength of the x rotation, while the frequency de- 
tuning, S = u — cjoi, gives the strength of the z rota- 
tion. As seen in Fig. [SJ^b), as the power increases, ef- 
fects of S errors due to hyperfine gradients decrease. At 
ujoi/27r = 0.355 GHz, the nutation frequency scales with 
voltage as dQu/dVi -- 27r x 70 MHz/mV. This scaling 
increases with dJx/de^ which grows as the 111 region 
is shrunk (^q -^ 0) to increase cjoi- At cJoi/27r = 1.98 
GHz, this scaling was measured to be ~ 27r x 5 GHz/mV, 
demonstrating a way to increase coupling to external 
voltages. 

On resonance in the rotating frame, the Hamilton takes 
the form Hrf = cos{^)ax -\-sm{^) ay, where ^ is the rel- 
ative phase of the carrier wave with respect to the first 
pulse incident on the qubit. Controlling phase relative to 
the initial pulse thus allows full two-axis qubit control. 
To test the qubit response, we prepare a |0) and drive 
a rotation on resonance for a time Tx, then apply a sec- 
ond pulse at relative phase ^ to drive a 37r/2 rotation 
in a time 37r/2co'R. Figure 4 shows data for <l> = 0°, 90°, 
and 180°, along with model curves using an optimized, 
though reasonable, value for hyperfine couplings as a fit 



parameter. 

Phase control was sufficient to implement a CPMG dy- 
namical decoupling sequence, where 7r-pulses are applied 
along the y axis in the rotating frame, partially decou- 
pling rotation errors [15 . Figure [s] shows resulting coher- 
ence time, T2, for CPMG sequences up to 64 7r-pulses, 
which gave T2 = 19 ± 2 /is. Values for T2 were extracted 
from Gaussian fits to Po(rD), where td is the total de- 
phasing time (see inset of Fig. 4). Between 2 and 16 
pulses, the scaling of coherence time with (even) pulse 
number, n^^, appears well described by the power-law, 
T2 = A(n^)^, where 7 = 0.84 ± 0.05. Within a classi- 



cal power-law noise model 



implies S{uj) 



with a /3 = 5±1. The inconsistency of this result with re- 
cent studies of electrical noise in the singlet-triplet qubit, 
where /3 ~ 0.7 [26 , may reflect first-order insensitivity of 
the resonant exchange qubit to detuning noise. However, 
a detailed model for dynamical decoupling that distin- 
guishes voltage noise from hyperfine noise has not been 
developed to date. Moreover, pulse sequences designed 
to decouple hyperfine noise for exchange-only qubits [27^ 
may also be adaptable to the resonant exchange qubit. 

For n-j^ > 16, T2{n-j^) falls below the steep power-law, 
and appears to saturate around 20 jas. The measured 
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FIG. 5: T2 for various orders of CPMG-n, were each se- 
quence contains n tt rotations about ^, as depicted in the 
lower inset. The upper inset depicts the detuning sequence 
for this experiment. We found that up to n = 16, the even 
number of pulses was well described by T2 = ^(717^)^, where 
7 = 0.84 zb 0.05. This translates to a power spectral density 



of 5'(cj) 



\ where/3 = 5=bl. 



Ti for a splitting co'oi/27r = 0.33 GHz was ~ 40 /is, and 
decreased monotonically with increasing cjqi, consistent 
with phonon-based relaxation, which suggests that Ti 
was not limiting T2 at cjoi/^tt = 0.2 GHz. Pulse er- 
rors are likely limiting T2 in this measurement, though 
extending coherence much longer will require extending 
Ti. 

In summary, we have introduced and demonstrated the 
operation of a new quantum-dot-based qubit that uses 
exchange for both the longitudinal and oscillatory trans- 
verse fields. A large exchange gap prevents state leakage, 
and the operating point is insensitive to first order to 
fluctuations in gate-controlled detuning. Two-axis con- 
trol and a large ratio (~ 10^) of coherence time to gate 
operation time were demonstrated. Implementation of a 
two-qubit gate [161 EH] is next experimental challenge. 
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Supplementary Information for The Resonant Exchange Qubit 
MEASUREMENT AND NORMALIZATION 



A single normalization procedure was used for all data in the main paper to convert the measured reflectometry 
signals into output probabilities. It is similar to the normalization procedure described in Ref. [1 , except that the 
normalization of the Ti decay does not depend on measuring the overlap between \Sr) and l^^). A measurement 
of this is impossible because the overlap |(S'r|6'^) p requires diabatic passage through the center of the 111 region, 
and in our current setup the large gap in the center region forces our state to adiabatically follow the lower branch. 
Normalization of Ti is done through a separate independent measurement of Ti for each figure, as described below. 

As in Ref. [1 , the measurements of a given parameter (e, a;, burst power, etc.) was repeated 2^^ or 2^^ times to 
obtain measurement statistics and then histogrammed, following the procedure in Ref. |2|. The resulting histogram 
is fit to a function of the form. 



n(Vrf) = 



V^ 



exp 



(Vri-Vin 



102^2 



2a2 



-6"^^/^^^^=^ exp 



V^ 



(Vrf-V^n 



111^2 



/ 






(N tm {1-P ) 

V27rcr2 Ti A^rf 



2a2 



exp 



tm V- 



,102 



Kf-v)2 



Ti Avrf 



2(j2 



(SI) 



where n{v^i) is the fraction of histogram events with outcomes v^f for a measurement in 102, v]^ is the refiected 



voltage corresponding to double occupancy in the right dot. 



,111 



is the refiected voltage corresponding to single 



charge occupancy in the all three dots, Av^f 



,102 



, P is the fraction of 102 outcomes in the data set, Ti is the 



relaxation time at at the measurement detuning £^, tm is the measurement time, and a is the standard deviation of 
the histogram peaks due to noise in the rf equipment and shot noise intrinsic to the rf sensor dot. 
The extracted parameters v]^'^ and v]l^ are then used to normalize the return probabilities P as 



P = 



(^rf) 



,111 



,102 



,,111 



(S2) 



where {v^i) is the average voltage for a particular parameter (e, uo^ burst power, etc.) over all repetitions of the 
measurement sequence. 



Equation (S2) converts v^f into a probability, but it does not account for relaxation during the measurement time 
Tm, where a 111 state relaxes to a 102 state. As described in Ref. [3^, P is related to the actual probability Pq through 



^M Jo 



dt P exp 



t 



= P 



P\ 



1 — exp 



Tm 
Ti 



(S3) 



By knowing tm and Ti , we can correct for measurement relaxation. A measurement of the relaxation time at e^ 
acquired for each section of data in the main paper, by fitting the average probability as a function of tm- 



IS 



MODEL OF THE EXCHANGE INTERACTIONS 



We find that the detuning dependence of the exchange interactions Ji and Jr are well described by the model of 

(S4) 



a , 



Ji = -2(^+^0) + 



\ 



texp 




a^ 



+ Y(^+^o)' 



a 



(e 



^OJ + A, 



iexp 



WteJ 



-{e 



(S5) 



where a is the lever arm between e and energy, Wt is a phenomenological suppression of the tunnel coupling with e, 
and ±£o is ^^^ detuning of the 111-102 and 111-201 charge transitions, or half the width of the 111 region. We can 
then write Jz = \ {Ji{^) + Jr{^)) as 



J. 



I 



OLS^ - A 



texp 




+ y(£+£oP 



\ 



iexp 






+ Y(^-^o)' 



(S6) 



At £ = 0, eq (S6) simplifies to 



K,, - Vl 



\ 



texp 



1 



n \ 2 



+ Y^(^™-K^)^ 



(S7) 



where we have replaced e^ with (F^ — V^)/2. Experimentally, we see that the width of the 111 region, 2£q, is linear in 
Kn — K^, with the same lever arm as the other gates. Equation (S7) is used in Fig. 2(d) of the main text to map the 
resonance as a function of V^a [4 • In Fig. 2(d), t = 16.9 /ieV, Wt was taken to be very large, such that the exponential 
was ignored, V^ was taken to be -4.05 mV on this plot [5 . 

The transverse exchange, Jx = ^ (^r(^) — Ji{^))^ can be written as 




1 



texp 




+ ^{e+e,f 



texp 



\Wteo 



+ ^{e-e,f 



(S8) 



MODEL USED IN FIG. 2(C) OF THE MAIN TEXT 



Reference [6] gives tfie separation between tfie lower branch of the qubit state, which they refer to as |A'), and the 
\Q) state as 



Ea'Q = -^[Jl + Jr+ ^Jf + J^2 - JiJ^ ) , 



(S9) 



which is the difference between the lowest two qubit eigenvalues of equation (S17). The separation between \Q) and 
|(5+), Eqqj^ is g* jiBQ^f In Fig- 2(c) we plot the intersection of these two curves, 5ext(^) = ^o + ^A'q(^)/%*A^, 
where Bq is an experimentally determined offset in the field due to remnant fields from ferromagnetic components in 
the cryostat. We find that an offset of Bq = —9.3 mT and an effective g-factor of ^* = —0.34 describe our data well 
in the center of 111. The tunnel coupling t was 16.9 /ieV, Eq = 3.7 mV, a = 40/ieV/mV, Wt = 3. 



MODEL AND POWER BROADENING IN FIG. 2(D) 



The model in Fig. 2(d) is a plot of eq. (jSTf, where S = {Vi - V^°) + {Vr - V^) + j{Vm - Vm^) was held constant, 
with an experimentally determined 7 = 3. In Fig. |Sl[b), the resonance cjqi was extracted along with its width in 
frequency space by fitting it to a Gaussian at each value of Vm- We find in Fig. [sTj^d) that the resonance width, 
shown in Fig. [sTj^c), is proportional to duuoi/dVm, which could suggest that the resonance widens with electrical 
noise. duJoi/dVm is also proportional to dJ^/ds, which sets the strength of the Rabi oscillation. From this, we cannot 
determine whether the resonance is broadened due to fluctuations in Vm^ or due to power dependent broadening from 
an increased dJx/de. 
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Figure SI: (a) The data in F ig. 2(d) without the theory overlay, (b) The resonance center value extracted from Gaussian 
fits, along with the plot of eq. ( |S7[ ). (c) The widths of the resonance extracted from the Gaussian fits, (d) The widths of 
the resonance plotted against the analytic derivative of the qubit splitting uj — Jz/h^ where JziVm) is given by eq. (S7). The 
analytic derivative was chosen to remove the noise from the numerical derivative of the data in panel (b) . 



MODEL USED IN FIG. 2(F) OF THE MAIN TEXT 

The model in Fig. 2(f) used seven of the eight spin states available to the three electron system to reproduce the 
spectroscopic data. 
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The states beyond the two levels of the qubit manifold are included to account for the possibility of nuclear 
mediated leakage from the qubit subspace < Q+s \ , Q^i ) , Q_k y > as well as accidental initialization into the two 

states < 0-1/ •> 1-1/ f ^^^^ ^^^^ ^^^ same total spin as the qubit states, Sz = 1/2, but opposite spin projection, 
rriz = —1/2, as shown in Fig. S2 In the regime that the device is operated in, the Zeeman splitting due to the applied 
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Figure S2: The full energy level spectrum of the three electron system. The dashed lines are states with opposite spin 
projection to those displayed in the main text. 



magnetic field is less then the electron temperature, preventing us from preferentially loading 0^ i y and 1+ 1 ) instead 
of 



0_ 



and 



. The simulations presented in Fig. 



electrons from the leads we can in fact prepare 



^0_i), l_i), Q_i ) f i^ subsequent simulations. 
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0^1 y and 



give some indication that by avoiding replenishing our 
, which allows us to ignore the higher energy states, 



The eighth level, \Q_3) = |^^^), is separated from all other levels by the external magnetic field as shown in 

Fig. |S2[ and is therefore ignored in order to speed up computation. The simulation described the qubit evolution in 
the presence of exchange interactions and Zeeman energy from longitudinal and transverse nuclei, which we account 
for with the following Hamiltonians written in the basis of 
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Here, g* ^ —0.34, as determined from Figs. 2(c,e). The magnetic field terms Bxi^ Byi^ and Bzi in equations ( |S19[ ), 
(S20), and (S18) are the magnetic fields along f, ^, z respectively in dot z, where i = 1 corresponds to the left. The 



exchange terms Ji and Jr in eq. (S17) are the e-dependent terms from eqs. (S4) and (S5) 



The model was created in the following way. At a given e and cj, nine random variables were drawn from a 
normal distribution to take the nine nuclear field components, B^^ By^ and Bz in each of the three dots. From there, 
eigenstates of the full Hamiltonian, H = Hj -\- Hbx -^ Hsy -\- Hbz^ were calculated, and an initial state density matrix 
was chosen as a mixture of 90% of the eigenstate with the largest overlap with 0,1), and 5% of the eigenstates with 



the largest overlaps with 



+ 4 



and 



Q I 1 ) to account for loading infidelity. The initial state was then time evolved 



according to the Liouville-von Neumann equation. 



ih 



dp 



[H,p], 



(S21) 



for 300 ns in the presence of an oscillatory s. The final density matrix was then transformed into eigenstates of 
only the exchange interactions and the external magnetic field, and the population of the lower qubit eigenstate was 
recorded [?]• This process is then repeated 25 times with new random values for all of the nuclear fields, and the 
average return probability is recorded in the model ^. 

In the model, the amplitude was a 0.225 mV oscillation in detuning, equivalent to a 0.45 mV oscillation in Vi, or 
—51 dBm. The standard deviation of nuclear gradients was 3.9 mT with a ^* = —0.34. The tunnel coupling t was 
16.9 /ieV, £0 = 3.7 mV, a = 40/ieV/mV, Wt = 3. 

The simulation used seven levels to enable us to check whether we loaded in states from the ruz = — 1/2 manifold 

\ 0_i y , 1-1 ) , Q-i ) (' Figure S3^a) shows a detailed view of the region near the Q_^3 ) anti-crossing from the 
simulation Fig. 2(f) in the main text. Panel (b) shows that same region with a mixture of both rriz = +1/2 and 
rriz = —1/2 states loaded into the initial state, approximating the density matrix in the case that we load into the 
other Zeeman manifold approximately half the time. The presence of these extra states causes a doubling of the 
transition lines near this anti-crossing, a phenomenon that we do not see reproduced in the data, shown in Fig.|S3[c). 
From this we determine that we are only loading in the rriz = +1/2 manifold. This is consistent with our initialization 
procedure which maintains isolation from the higher temperature leads. 
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Figure S3: A detailed view of the spectroscopy described in the main text in Figs. 2(e,f), near the anti-crossing between 
the lower qubit branch and Q^s y at £ ~ —3 mV. (a) A simulation of the evolution of our three electron system near 

the Q^3 ) anti-crossing, where the initial state was contained in the niz = +1/2 manifold s O^i ) , l^i ) , Q+i ) \- The 



slanting bright lines indicate transitions mediated by n-photon of frequency u between the lower qubit branch and Q^s 

while the vertical bright line is the anti-crossing itself. Transitions to the right of the anti-crossing represent the stimulated 
emission of photons, while transitions to the left indicate the absorption of photons [see Fig. 2(a) of the main text and 
Fig. |S2| . (b) A simulation where the initial state was contained states from both the ruz = +1/2 and ruz = — 1/2 manifolds 

i I 1 ) , 1 I 1 ) , Q I 1 ),0_i),l_i), \Q_i ) >. Here, the slanting bright lines are doubled, indicating an additional set of 

L|"'"2/ |"'"2/ I "'"2/ I 2/ I 2/ I 2/J 

photon mediated transitions which are slightly offset in detuning 9 . This implies that the presence of iriz — —1/2 states in 

the initial state should appear in the data as a double of the spectra around the Q^s ) anti-crossing, (c) A detailed look at 

the data corresponding to the same range of the simulation. The anti-crossing is shifted slightly from e ^ —3 mV due to small 
discrepancies between the exchange model and the data at large exchanges. The width of the photon lines in (c) is less than the 
spacing between the doubled lines in (b), leading us to conclude that rUz = —1/2 states were not loaded in significant amounts. 
A more detailed measurement which includes a spin selective readout of the third electron 10 would be necessary to confirm 
this result conclusively. 



MODEL USED IN FIG. 3 INSETS OF THE MAIN TEXT 



Having shown in the previous section that we are only in the space of O^i V l+i ), and 



^1 ) states, we now 

restrict ourselves to only the appropriate 3x3 subregions of Hamiltonians (S17) and (S18), which is valid when we 
are far away from the \Q-\-) anti-crossings as we are in Fig. 3. We express these as 
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(S23) 



with Hbx = 0, and HBy = 0. We also revert to the simpler notation of |0), |1), and \Q). For the insets, we create 
an initialization state that is 90% |0), and 10% |1) and |Q), and time evolve it using equation ( |S21 ) for 100 ns in the 
presence of an oscillatory e and static longitudinal nuclear field gradients ABi and ABr. The |0)(0| term is extracted 
at multiple times during the evolution, and recorded. This process is repeated with a new set of ABi and ABr drawn 
from a random normal distribution, and averaged with the previous set. Unlike the simulation in Fig. 2(f), the nuclei 
are presumed to be static over the course of a single column, which provides the flickering effect seen in the model 
and the data. The data was acquired more rapidly in this set, justifying the modification to the model. 

In the model in Fig. 3(a), the amplitude was a 0.45 mV oscillation in detuning, equivalent to a 0.89 mV oscillation 
in l^, or —45 dBm. The standard deviation of nuclear gradients was 3.9 mT with a ^* = —0.34. The tunnel coupling 
t was 16.9 /ieV, £q = 3.7 mV, a = 40/ieV/mV, Wt = 3. 



12 
MODEL USED IN FIG. 4 OF THE MAIN TEXT 

The theory curves in Fig. 4 are based on a similar model to the insets in Fig. 3, with the complication that we 
perform a second time evolution which takes the final state of the first evolution as its input. In the second time 
evolution, the e oscillation at a phase ^ with respect to the oscillation in the first time evolution. The preparation 
was 90% |0) and 10% |1) and IQ), and both longitudinal as well as transverse nuclear fiuctuations were incorporated 
in a quasi- static manner, as they were in the inset of Fig. 3. 

In the model curves in Fig. 4(b), the amplitude was a 0.30 mV oscillation in detuning, equivalent to a 0.60 mV 
oscillation in Vi^ or r^ —55 dBm. The standard deviation of nuclear gradients was 5.2 mT with a ^* = —0.34. The 
tunnel coupling t was 12.4 /ieV, £q = 3.7 mV, a = 40/ieV/mV, Wt = 3. The somewhat larger standard deviation 
of nuclear gradients in this data compared to the data in Fig. 3 is not well understood, but since the nuclear field 
fiuctuates slowly on the order of minutes, we expect to observe a range of standard deviations among individual data 
sets. 



[1] J. Medford, J. Beil, J. M. Taylor, S.B. Bartlett, A. C. Doherty E. I. Rashba, D. P. DiVincenzo, H. Lu, A. C. Gossard, 

and C. M. Marcus, arXiv preprint arXiv: 1302. 1933 cond- mat. mes- hall] (2013). 
[2] C. Barthel, D. J. Reilly C. M. Marcus, M. P. Hanson, and A. C. Gossard, Phys. Rev. Lett. 103, 160503 (2009). 
[3] A. C. Johnson, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 

Nature 435, 925 (2005). 
[4] In Fig. 2(d), the tunnel coupling was simplified back to just t, with no exponential dependence. The physical reason for 

this is not well understood, but the model clearly fits the data if this is true, and drastically disagrees if the exponential 

dependence is included. 
[5] An overall background to Vm of —452 mV was already removed, leaving just the amount that was pulsed using the AWG. 
[6] E. A. Laird, J. M. Taylor, D. P. DiVincenzo, C. M. Marcus, M. P. Hanson, and A. C. Gossard, Phys. Rev. B 82, 075403 

(2010). 
[7] The eigenstates of the Hamiltonian with the nuclei were used for the initialization because the slow ramp in to this state 

allows the system to adiabatically enter the anti-crossing with Qs ) . The pulse to measurement position is a much more 

I 2/ + 

rapid, meaning that the qubit left that region diabatically The difference in eigenstates allows us to capture this effect 

reasonably well without considering the full time dynamics of the problem. 
[8] In order to speed up computation, the repetitions were truncated at 5 samples if the average return probability was > 86%. 

This saw a drastic improvement in speed with little qualitative change to the plots. 
[9] Panel (b) was sampled more sparsely because it was a test run of the larger initialization space. Since it did not accurately 

represent the data, it did not receive a follow up calculation at higher sampling, as panel (a) did. 
[10] J. M. Elzerman, R. Hanson, L. H. Willems van Bereven, B. Witkamp, L. M. K. Vandersypen, and L. P. Kouwenhoven, 

Nature 430, 431 (2004). 



